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The SU(2)l ® SU(2)r, symmetric Yukawa model with mirror-fermions in the limit where the mirror-fermion is 
decoupled is studied both analytically and numerically. The bare scalar self-coupling A is fixed at zero and infinity. 
The phase structure is explored and the relevant phase transition is found to be consistent with a second order 
one. The fermionic mass spectrum close to that transition is discussed and a first non-perturbative estimate of the 
influence of fermions on the upper and lower bounds on the renormalized scalar self-coupling is given. Numerical 
results are confronted with perturbative predictions. 



1. INTRODUCTION 

This contribution represents the talks "Phase 
Structure of a Chiral SU(2) Yukawa Model", 
"Mass Spectrum and Bounds on the Couplings in 
Yukawa Models with Mirror-Fermions" delivered 
at the international conference on lattice field the- 
ory "Lattice 92", Amsterdam, 15-19 September 
1992, by Lee Lin and Hartmut Wittig, respec- 
tively. 

The minimal Standard Model (SM) has been 
very successful. So far all the experimental data 
are in good agreement with its perturbative pre- 
dictions, but there are still too many free pa- 
rameters (18, if all neutrinos are massless) in the 
SM, hence physicists believe that it is at most a 
good low energy effective theory, and new physics 
which cannot be described by the SM will have 
noticeable effects at a higher energy scale. Here, 
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we concentrate on two free parameters in the SM: 
the top quark mass and the mass of the Higgs 
particle which is the remnant scalar particle after 
spontaneous symmetry breaking. We would like 
to investigate whether these two parameters are 
completely free, or whether there are some limi- 
tations on their values. Hopefully, one can also 
get some hints about the new physics by studying 
this Higgs-fermion sector in the SM. 

Since the presently quoted upper bound on the 
top quark mass of 200 GeV is a one-loop per- 
turbative result, nonperturbative (i.e. lattice) 
studies of this issue can be very helpful. In order 
to deal with the problem of fermion doublers in 
lattice formulations, we take the mirror-fermion 
approach 

Our calculations were done in the SU(2) version 
of the Yukawa model with explicit mirror pairs of 
fermion doublet fields. The lattice action is a sum 

of the 0(4) (= SU(2)l«)SU(2)r) symmetric pure 
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scalar part S^p and the fermionic part Sx^: 

S = Sip + . (1) 

ipx is the 2 (g) 2 matrix scalar field, and '^x = 
a; J Xx) stands for the mirror pair of fermion dou- 
blet fields where "0 is the fermion doublet and x 
the mirror- fermion doublet. In the usual normal- 
ization conventions for numerical simulations we 
have 



±4 . 

- K ^ Tr {vt+[.'Px) \ ' 

u=±l ^ 



(2) 

X ^ 

±4 

[(Xx+A'/':!;) - (Xa;^'^;) + {i^x+i,Xx) - i^xXx)] 
+ [{ll^BxVt'^Lx) + {lpLxfxll^B,x)\ 

+ [{XRxfxXLx) + (XLa;<^^Xi?:E)] ^ • (3) 



Here K is the fermion hopping parameter, r the 
Wilson-parameter, which will always be fixed to 
r — 1, the indices L, R denote the chiral compo- 
nents of fermion fields, and Tr acts on the SU(2) 
indices only. In this normalization the fermion 
mirror- fermion mixing mass is /z^^; = 1 ^ %rK . 
The lattice spacing a is set to unity. 

At Gx = 0, the action has the Golterman- 
Petcher shift-symmetry such that all higher ver- 
tex functions containing the X'fisld vanish iden- 
tically d,^. In this hmit, the X'X and X'i' com- 
ponents of the two-point fermion vertex function 
(the inverse fermion propagator) {p) are equal 
to the corresponding components of the free in- 
verse propagator By setting K — 1/8, one can 
easily show that in the broken phase, there is no 
mixing between the fermion and mirror-fermion 
and the mirror-fermion completely decouples like 
a right-handed neutrino. This = K ~ 



1/8 combination is sometimes called the mirror- 
fermion decoupling limit. This leaves us only one 
mass parameter n to tune when we try to ap- 
proach the continuum limit and hence saves a lot 
of CPU time. Numerical simulations can also be 
guided by analytic results derived in this limit 
such that fluctuations can hopefully be reduced. 
The Hybrid Monte Carlo algorithm |^ is used in 
our simulations. We therefore need to double the 
flavour of the fermion spectrum to guarantee a 
positive definite fermion matrix determinant. 

One of the main goals in the analysis of the 
model is the determination of the range of renor- 
malized quartic couplings gn that can be realized 
for a certain value of the cut-off as a function of 
the renormalized Yukawa coupling Grtp . This so- 
called "Allowed Region" is bounded by the (up- 
per) triviality bound on (saturated at A = oo) 
and the (lower) vacuum stability bound, which 
is defined at A = 0. Hence for the bare quartic 
coupling A the values A = cxd, are crucial in the 
numerical simulations. 

Reflection positivity, which is required for any 
Euclidean quantum field theory respecting uni- 
tarity, can be proven to hold in a wide range of 
the parameter space of our model . In particu- 
lar, reflexion positivity only holds for k > 0. Al- 
though the criterion of reflexion positivity is only 
a sufficient but not a necessary condition for the 
reconstruction of the theory in Minkowski space, 
we would like to stay where unitarity is guaran- 
teed in order to be on the safe side. So we always 
study the "Allowed Region" for the renormalized 
couplings in the subspace where k is non-negative. 

The /3-functions have been calculated up to one 
loop on the lattice and two loops in the contin- 
uum Q. There are basically two possibilities 
for the cut-off dependence of the "Allowed Re- 
gion" : either this region shrinks to the origin as 
the cut-off grows or it expands and eventually fills 
the whole space of renormalized couplings at in- 
finite cut-off. The first possibility is due to the 
fact that there is only one fixed point at zero 
couplings and that fixed point is infrared stable. 
One-loop /3-functions are behaving qualitatively 
like this. In the second scenario, there is an ul- 
traviolet stable fixed point at nonzero couplings 
which changes the dependence of the "Allowed 
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Region" on the cut-ofF qualitatively. Assuming 
that the /3-functions are behaving qualitatively 
like at one-loop level, then our model is trivial in 
the continuum limit, and the upper bound on the 
renormalized scalar self-coupling gn will be given 
by setting A = oo while the lower bound comes 
from the A ^ limit . 

2. PHASE STRUCTURE 

In order to know where and how the cut-off 
can be removed, we need to explore the phase 
structure. If there is a first order phase transition 
somewhere in the bare parameter space, we would 
like to know how it is going to affect the "Allowed 
Region" . 

We basically use the magnetization {\(f)\) and 
staggered magnetization as the order param- 
eters to distinguish different phases where ( ) is 
the statistical average, and 

|</)| ^ , 101 ^ , 

X 

X 

where L and T are lattice sizes along the spa- 
tial and time directions. The symmetry bro- 
ken (FM) phase has nonzero (|(/)|) and vanishing 
the symmetric (PM) phase has zero (101) 
and (101) while (|0|) = 0, (|0|) 7^ in the anti- 
ferromagnetic (AFM) phase. In the PM phase, 
the mirror- fermion and fermion are degenerate 
and hence it is not physical. The physically rel- 
evant phase is the FM phase where the split- 
ting between fermion and mirror-fermion masses 
is possible due to spontaneous symmetry breaking 
and by tuning the bare couplings appropriately. 
The AFM phase is like the FM phase except that 
the staggered scalar field is now playing the role 
of 0. 

The phase structure can be investigated ana- 
lytically using various expansions in several lim- 
its 10]. We find that at zero or infinite Yukawa 
couphngs, or at if = 00, the system goes to a 
purely scalar four-component 0^ theory plus free 



fermions. Therefore, the system exists in the FM, 
PM and AFM phases with the transitions be- 
tween them being second order and Gaussian. At 
X = 0, which is the limit where the bare fermion 
mass is infinite, a first order phase transition line 
was found in the U(l) version at • > and 
finite A |^ . This first order transition is due to the 
singularity of a log-term in the effective action. It 
becomes weaker and weaker and eventually van- 
ishes at A = 00. In the SU(2) version, everything 
is qualitatively the same. Therefore, this first or- 
der transition is also present when G^ ■ G-^ > 0, 
A = finite. We did not further investigate this is- 
sue. At A = cx), when Yukawa couplings are weak 
or strong, one can do expansions up to the next- 
leading order in \G\ or 1/|G| plus the small-K 
expansion and find out that Kc, the value where 
the transition between the FM and PM phases 
happens, will go down as \G\ or 1/\G\ increase. 
Furthermore the transition remains a second or- 
der one on which the renormalized scalar mass 
vanishes. 

We also need to look for the fermionic critical 
plane on which the renormalized fermionic mix- 
ing mass ^.R vanishes. When = at any 
A, this happens at K ^ \/S and any G^ value. 
When both G's are nonzero, we can use one-loop 
bare perturbation theory to estimate its position 
at weak K and G's. 

The correct continuum limit should be taken 
by approaching the critical line where both scalar 
and fermion masses are zero from within the FM 
phase while at the same time the mixing mass is 
fixed at zero and all mass ratios are kept constant. 

The phase structure at intermediate values of 
the Yukawa couplings must be explored numeri- 
cally. That was always done on the 4^ • 8 lattice. 
The strategy is to fix A at infinity or zero (actually 
10~^), which are the relevant values to studying 
upper and lower limits on t^^; respectively. For the 
remaining parameters, since we chose to decouple 
the mirror-fermion like a right-handed neutrino in 
the SU(2) model, we always set G^ = 0, X = 1/8, 
and then studied the phase diagram in the (k, G^) 
plane. 

At A = cx), the phase structure was explored 
numerically up to G^p around 2.0. The result 
was shown in Fig.l in ref. Q|. Everything ap- 
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peared to be qualitatively the same as the phase 
structure of the U(l) model at infinite A. At 
> 1.5, K, = —0.5, a new phase with nonva- 
nishing magnetization and staggered magnetiza- 
tion was found. We call it the ferri- magnetic (FI) 
phase. The phase transition between the FM and 
PM phases was found to be consistent with a sec- 
ond order one because of the smooth behaviour of 
(|(/)|) across the transition. This was also true for 
other transitions in the phase diagram. Based on 
our analytic analyses, we know that as we move 
to larger and larger values, eventually the FI 
phase will disappear and the system again exists 
only in the FM, PM and AFM phases. We there- 
fore did not spend CPU time to investigate fur- 
ther. (This has been numerically confirmed in our 
U(l) model §.) 

The phase structure at very small A and large 
negative values of k was actually not thoroughly 
explored yet. In the U(l) model we only made 
sure that the FM-PM transition was consistent 
with a second order one and then went on to study 
the lower bound on gj^ We now spent more 
computer time to look into the phase structure of 
the SU(2) model at A = 10"^. The transitions 
from the FM to PM and from the PM to AFM 
phases were found to be consistent with a second 
order one at weak G^. Again the magnetization 
behaves smoothly across the PM to FM transi- 
tion. Plotting histograms of the magnetization 
we found that there is no evidence at all for a two- 
state signal. The two transition lines bend down 
as G^ increases, and come quite close to each 
other at G^ = 0.75, k, = —0.17. At an even more 
negative K-value: —0.19 and at G^ = 0.6,0.8, 
we observed that the system seemed to "tunnel" 
from the FI phase to the AFM phase. At this 
stage, we cannot decide whether it is a real tun- 
neling or whether the system has not equilibrated 
yet. Even if we find a tunneling event (or a hys- 
teresis loop), it still might be due to a very long 
auto-correlation length in the vicinity of a critical 
point and therefore is still not the final word. We 
think it is better to measure the shape of the ef- 
fective potential to see if a double-well structure 
developes. If there is a first order phase tran- 
sition, we suspect that it might be the continu- 
ation of the first order transition we discovered 



at K = 0. On the other hand, a leading order 
large- calculation of the effective potential of 
our SU(2) model shows that the usual PM, FM, 
AFM and FI phases exist at small A, and that all 
transitions are second order due to the absence of 
a quartic term in the effective potential in leading 
order. This issue will be further investigated in 
the future. 

Most importantly, the physically relevant FM- 
PM phase transition at A = 10~^ is consistent 
with a second order one, and we can define the 
continuum limit by approaching it from the FM 
phase. Therefore, studies of the lower limit on gn 
should not be affected by a possibly existing first 
order transition. 

3. MASSES AND COUPLINGS 

We now describe the Monte Carlo simulations 
of the model in the broken (FM) phase. Besides 
the massive component (Jx of the scalar field <Px 
three massless Goldstone bosons TTja,, j = 1,2,3 
appear which cause the strong finite-size effects 
encountered in the simulations. 

As noted before, the model was simulated us- 
ing the Hybrid Monte Carlo algorithm. For the 
fermions, periodic spatial boundary conditions 
but antiperiodic boundary conditions in the time 
direction were chosen. This fixes the minimum 
lattice momentum for fermionic quantities at 
Pmin = (0,0, 0,7r/T) where T is the time extent 
of the lattice. Exploiting the shift symmetry, the 
decoupling of the mirror-fermion was ensured by 
setting — 0. In order to exclude mixing be- 
tween ip- and x-states we chose a slightly different 
value than K — which was originally sug- 
gested. Namely, by setting /ip^j,, = 0, where 
l^p = {iJ-ijix + KrfP') /2K , the fermionic hopping 
parameter K turns out slightly greater than 1/8, 
but will eventually approach this value once the 
time extent of the lattice goes to infinity. This 
particular choice corresponds to exact decoupling 
in the continuum limit whilst ensuring a smooth 
behaviour of the propagator on a finite lattice 
near if = 1/8 (see also ref. 0] ) . 

With the parameters G^ and K fixed by the 
decoupling condition our further strategy was as 
follows. For A = oo, G^ was varied from 0.3, 0.6 
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to 1.0. At each value of the scalar hopping 
parameter k was tuned in order to achieve < 
1. Choosing lattice sizes of • T of 4^ • 8, 6^ • 12 
and 8'^ • 16 it was hoped that on the largest lattice 
niRa could be brought down to about mua — 0.5 
in lattice units. 

For A = 10^^, Gjf, was fixed at 0.3 and again 
K was tuned. A second set of data, however, was 
generated by setting k — and tuning in or- 
der to determine the maximum value of the renor- 
malized Yukawa coupling Gr^ for non-negative 
values of k, i.e. in the region where reflexion pos- 
itivity can still be proven. 

The results for the scalar mass show that mucr 
is decreasing as one approaches the critical line 
from above and then starts rising again. This 
rise close to the phase transition is explained as 
a strong finite-size effect. The minimum of 
increases with increasing G^ [Q. For A = oo, 
= 0.6, it is necessary to have lattices as large 
as 8'^ • 16 in order to obtain scalar masses smaller 
than 1 in lattice units. 

At A = lO^*' the tuning of m/j^ is easier in the 
sense that one can achieve smaller values than at 
A = oo, however the rise of the minimum of 
as Gjp is increased, remains. All that illustrates 
that there are strong finite-size effects which get 
stronger as the bare couplings A and/or are in- 
creased. Nevertheless, our experience has shown 
that for K > a lattice size of 8'^ • 16 is sufficient 
to have control on the finite volume effects. 

For the fermion masses one can make the fol- 
lowing statements: firstly, setting G^^ = ensured 
indeed that the mass of the mirror-fermion fi^^ 
was always zero within errors. Secondly, the con- 
dition /ip„i„ = guaranteed that the fermionic 
mixing mass fj^R was zero as well. Fig. 1 shows 
the fermion mass fj,R^, (full squares) and the low- 
est doubler states (open squares) and fif/^ 
(open circles) as a function of k. It is seen that 
the lowest fermion doublers have masses of 1.5 in 
lattice units. Furthermore, it is evident from the 
figure that only at the lowest value of k there is 
a clear separation of fiR^ and the doubler states. 
Since this K-value coincided with the minimum 
of the scalar mass {niRa- = 0.75(3) at k = 0.09, 
A = 10~^ G^ = 0.3 on 6^ • 12), it is precisely 
this data point which was subsequently included 



in the plot of the Allowed Region. To summarize, 
we note that in general fj,R^ decreases as k and/or 
G^ is decreased, and that the lowest doubler state 
of the mirror-fermion has a mass of about 1.5, ir- 
respective of the actual choice of the tunable bare 
parameters. 

Using the relation = Gr^pvr, where vr is 
the VEV of the scalar field, it is instructive to 
plot Gr^ versus G^ in order to check whether 
the linear rise of Gr^ reported for the symmet- 
ric phase 1^ prevails. This is done in Fig. 2. 
It is seen that there is still a linear behaviour 
which, however, is much flatter as in the sym- 
metric phase. The data obtained at A = 10^^ are 
slightly larger than those at A = oo. The maxi- 
mum value of Gr^ for k > was found at k = 0, 
Gv; = 0.63 where G'j^f = 3.5 ± 0, 4. This has 
to be confronted with the tree unitarity bound 
which yields Gr^ ~ 2.5 for Nf = 2. Hence the 
maximum value for Gr^ is not very much above 
this bound, suggesting that for k > the coupling 
Gr^ cannot grow indefinitely large. 

Fig. 3 shows our data for the Allowed Region, 
together with the perturbative expectation com- 
ing from the numerical integration of the one- 
loop /3-functions for cut-off's corresponding to 
mRa = 0.75 (dotted curve) and mRc = 1 (solid 
curve), respectively. All data were obtained on 
lattices of size 6'^ • 12 and 8'^ • 16. The point with 
the largest couplings was on 6^ • 12 at A = 10~^, 
G^ — 0.63, n — 0. Despite the large error 
bars the data are in remarkable agreement with 
the perturbative estimates even at values of Gr^ 
around or even above the tree unitarity limit. 
Therefore, the data support the view that the per- 
turbative behaviour of the bounds on the scalar 
coupling observed in pure (j)^ theory persists when 
fermions are included. 

Finally, Fig. 4 shows the mass ratio fJ-R^p/mRa 
as a function of Gt?,^, together with the pertur- 
bative prediction from the integration of the one- 
loop /3-functions for the upper and lower bound, 
respectively. It is seen that for larger values of 
Gfl0 (where there is only a weak cut-off depen- 
dence) the data are clustered around 0.7, which is 
close to the fixed point for the corresponding ra- 
tio of couplings, V^GRjp/y/gR — 0.759 . . . that is 
reached as the cut-off scale becomes infinite. The 
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Figure 1. The fermion mass fiji^ and the low- 
est doubler states plotted versus k at A = 10~®, 
= 0.3 on 6H2. 



Figure 3. The Allowed Region in the {gn, Gj^^)- 
plane together with the integration of the /3- 
functions for scalar mass m.Ra- = 1 (solid line) 
and miia- =0.75 (dotted line). 



Figure 2. The rcnormalizcd Yukawa coupling 
GriI, versus G^ for A = 00 (open circles) and 
A = (fuU squares) on 6^ • 12. 



Figure 4. The mass ratio fiR^/rriRcr versus Gr^ 
in comparison with 1-loop pcrturbative estimates 
for niRa = 0.75 (dotted curve), m/jo- = 1 (full 
curve) and niRc, = 1.25 (dashed curve). 
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fixed point is indicated by tlie iiorizontal line in 
the figure. For smaller Gr^ the data obtained 
for A = oo deviate from the fixed point in a fash- 
ion that is well described by the expected cut-off 
dependence of the upper bound for finite cut-off. 

4. CONCLUSIONS 

According to our results, the Higgs- Yukawa 
sector of the Standard Model can indeed be 
studied non-perturbatively following the mirror- 
fermion approach. The phase diagram exhibits 
a rich structure with the familiar PM, FM, 
AFM and FI phases encountered in other lattice 
Yukawa models as well. For both A = oo and A ~ 
the physically relevant phase transition from 
PM to FM is second order and hence permits tak- 
ing the continuum limit along this transition line. 
Exploiting the Golterman-Petcher shift symme- 
try, numerical simulations in the decoupling limit 
in the broken phase arc feasible and greatly facil- 
itated since only n remains to be tuned once 
is fixed. The resulting fermionic mass spectrum 
consists of a fcrmion of mass /i/j^, < 1, whereas 
Hn-^ and the mixing mass jin vanish identically. 
The doublers for both fermion and mirror-fermion 
receive masses of at least 1.5 in lattice units. 

The most severe limitations arise from finite- 
size effects on the scalar mass niRa in the broken 
phase. However, for an exploratory study such 
as this a lattice size of 8"^ • 16 appears sufficiently 
large in order to control finite-size effects. 

Restricting the physical analysis to data points 
with K > (i.e. data which always will per- 
mit a particle interpretation in Minkowski space- 
time), it is evident that the results for the cou- 
plings ga and Gr^ agree with perturbative es- 
timates. In particular the data for the Allowed 
Region suggest that the perturbative behaviour 
of the bounds on encountered in pure the- 
ory prevails once the effects of heavy fermions are 
taken into account. 

During the preparation of this paper, we have 
accumulated about 3500 trajectories for a run at 
A = 10-6, G^ = 0.25, = 0.0 and k = 0.099. 
It shows that with niRa = 0.624(56), jiR^ = 
0.413(12), we obtain qr = 16± 5, Gr^ = 1.56(7) 
and . = 2.43(22). Readers can see that the 



central value is slightly above the dotted one- 
loop curve we show in Fig. 3. Since the scalar 
mass in that data point turns out to be higher 
than niRcy = 0.75 for the dotted curve, our find- 
ings still support the perturbative scenario for the 
boimds on g^. Data with better statistics will be 
published in a forthcoming paper. 
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